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PAINTING NEW LINES: 

MAXIMIZING COLOR DIFFERENCE IN METRO MAPS 

SIMONE GRIFFIOEN AND ARTHEMY KISELEV 


Abstract. Each metro line usually has its own color on the map. For obvious rea¬ 
sons, these colors should be maximally different. Suppose a new metro line is built. 
We now explain a strategy to choose the color or several colors for the new lines, 
making them as different as possible from both the old ones and each other. This is 
illustrated by using the Moscow metro map. 


1. Introduction 

• “Have you seen the new terminal on the greenl" 

• “The red is overcrowded every morning at this time.” 

• “I live at the very end of orange, still OK with me: the escalator brings me right 
into the office.” 

• “I have to change to purple here. See you tomorrow!” 

Sure, it is all about metro. The subway. The metropolitan underground railway 
system. Each line has its own color and each color is the name of a line. A problem 
arises when another line is built: the old lines must keep their colors whereas the new 
line should have a color which looks maximally different from the old ones. It gets even 
more complicated if there are two or more new lines in a subway. Not only the new 
colors must be unlike the old ones, but also they must differ from each other as much 
as possible. 

In this paper we discuss a mathematical strategy to choose these colors. By the 
way, what new colors would you suggest for the real metro map in Figure 1? For our 
purposes this map is ideal: not only there are 14 colors in it, but also there will be 
several new lines built in that subway in the near future (according to [2]). Let us 
take this map to illustrate our strategy — yet the reader can apply the techniques to 
whathever other map he or she likes. 

2. Measuring the color dieeerence 

To maximize color differences - here, in metro maps - we need a way to measure them 
and argue about colors objectively. For this we describe them as tuples of numbers 
by using some coordinate system (or color space). For examble, the RGB color space 
is well-known: in it every color is encoded by its red, green, and blue components. 
Unfortunately, in the RGB space the Euclidean distance is not proportional to the 
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Figure 1. It is likely that you have already recognised the Moscow metro^ here. 

visible distance that people perceive in between colors. That is, the RGB space is not 
perceptually uniform — which we would like. 

In 1931, the Commission Internationale de I’EcIairage (CIE) did measurements on vi¬ 
sion of hundreds of humans. From those data, several color spaces were created, to make 
communication about colors easy and unambiguous. In 1976 the CIE recommended the 
use of two approximately perceptually uniform color spaces and color-difference formu¬ 
las. These spaces have become known by their officially recommended abbreviations 
CIELAB and CIELUV. In the former, a color is described by three coordinates: L 
(lightness), a (red-green scale) and b (yellow-blue scale).^ We will use this space to 
describe the colors of the Moscow metro map in terms of numbers: their CIELAB co¬ 
ordinates can be found in Table 1; Figure 2 shows a stereographical plot of the colors 
of the Moscow metro in the CIELAB space. 

The CIELAB space itself is not bounded; however, based on what one wants to 
measure, only part of it matters in practice (for example, only the colors that a printer 
can print or only the colors the human eye can see). Such a subset of colors is called 
a color gamut. The gamut of Adobe 98 in the CIELAB space can be seen in Figure 3. 
To keep the problem simple, we aproximate this gamut by a polyhedron. Because the 

^Opened in 1935 with 11 km of lines and 13 stations, the Moscow metro has 196 stations nowadays 
and the total route length of its lines is at least 327,5 km [4]. It is one of the world’s most heavily used 
metro systems. 

^The L-scale runs from 0 to 100; the parameters a and b typically have values in the range -150 
to 150. In this article we will round the coordinates of the colors to integers. The human eye can 
discriminate up to ten million colors, so it can distinguish the differences even finer than integer steps. 
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Figure 2. A stereographical image of the colors of the Moscow metro 
map at their location in the CIELAB space. The L-axis is the vertical 
axis, the a-axis runs from south-west to north-east and the 6-axis from 
south-east to north-west. {To see a 3-dimensional picture, view the left 
picture with your left eye and the right picture with your right eye at the 
same time.) 


CIELAB space was supposed to be almost perceptually uniform, the CIE 1976 color 
difference formula (CIE76) is given by the Euclidean distance 

AE,, = V(Li-L2)2 + (ai-02)2+ (61-62)2. 

In fact, the CIELAB space is not really uniform. In particular, at high values of a and 
6, the far too simple CIE76 formula values color differences too stronly compared to 
the experimental results on color perception. This is why several new color-difference 
formulas were proposed, of which the newest is the CIEDE2000 formula^ (or AEqo). 

3. Mathematical formulation of the problem 

Given the CIELAB space equipped with the Euclidean metric - or even better, the 
CIEDE2000 color difference formula - and the color gamut T (which we chose to be a 
polyhedron), we consider a set of existing colors P = {pi,... ,pk} within this gamut, 
represented by their CIELAB coordinates (e.g., see Table 1). If we want to add only 
one new line, this leads to the following maximin problem: 

Problem 1 (one new color). How can we find an x eT such that 

min lla; — Pi\\ —>■ max ? 

PiSP x&r 


^In the CIEDE2000 distance some weighting coefficients have to be chosen. We use the ratio 
= I : 1 : 1. This means that we consider the difference in chroma c = ^/aP+TP and 
hue h =arctan(6/a) to be more important than the difference in lightness L. For more information 
about this color difference formula, see [3]. 
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Figure 3. The gamut of Adobe 98 in the CIELAB space. We approx¬ 
imate this gamut by a polyhedron by taking the cornerpoints of this 
polyhedron close to the points A,... H in the hgure. 


Line 

Color 

L 

a 

b 

min AEqq 

1 

• 

Red 

52 

74 

53 

29,1 

2 

• 

Ocean green 

56 

-45 

26 

21,7 

3 

• 

Cobalt 

35 

7 

-43 

19,4 

4 

• 

Sky blue 

61 

-16 

-42 

6,3 

5 

• 

Olive brown 

41 

6 

27 

20,4 

6 

• 

Peach 

76 

24 

67 

14,5 

7 

• 

Pinkish purple 

43 

64 

-24 

28,3 

8 

• 

Light mustard 

86 

4 

85 

14,5 

9 

• 

Light grey 

71 

0 

-2 

9,4 

10 

• 

Greenish yellow 

80 

-26 

68 

17,8 

11 

• 

Turquoise blue 

56 

-21 

-29 

6,3 

Water 


Pale blue 

93 

-8 

-9 

11,3 

Background 

o 

White 

100 

0 

0 

9,4 

Text 

• 

Black 

13 

2 

0 

20,4 


Table 1. Coordinates of colors on the Moscow metro map. The names 
of the colors come from [6]. The average minimal distance is 16,3. Note 
that some of the minimal distnces are the same; this happens when the 
two points are the closest to each other. 


That is, the minimal distance from the new color x to the existing ones must be as large 
as possible, provided that the color x stays inside or on the boundary of the gamut T. 

If we want to add several new lines, each having its own new color, the maximin 
problem becomes this. 
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Problem 2 (m new colors). How can we find a collection ofm points X = {xi,..., Xm} 
such that 

min {\\xi - Xj\\,\\xi - pk\\} ^ max7 
Hi 

Pk&P 

That is, we maximize both the minimal distance in between the new colors xi,, Xm 
and their minimal distance to the old colors. Alternatively, we could add the new colors 
Ti, ... ,Xm one by one, using a method for solving Problem 1. However this does not 
always yield the optimal result.^ 

4. The Voronoi diagram method 

Let us consider Problem 1. Because the CIELAB space is continuous (at least, in 
principle), it is impossible to “try all points” in the gamut and hnd the one with the 
maximal minimal distance. Therefore it would be much easier to have only a hnite 
number of candidate points to check. This can be achieved by using Voronoi diagrams. 

Definition 1. Let P = {pi,... ,pk} C be a finite set {e.g., of old colors) such that 
Pi 7 ^ Pj for i j ■ The Voronoi polyhedron of a point Pi is the region 

V(pi) := {t G : ||a; — Pi\\ < ||a; — pj\\ for all j 7 ^ i}. 

The set V(P) = {V(pi),... ,V{pk)} is the {3-dimensional) Voronoi diagram generated 
by the set of points P. 

It is easily shown that the sought-for point x maximizing the minimal distance in 
Problem 1 belongs to the intersection of 

either: three Voronoi planes, that is, three faces of Voronoi polyhedra, 
or: two Voronoi planes and a boundary plane of the gamut, 
or: a Voronoi plane and two boundary planes, 
or: three boundary planes. 

This means that as soon as these candidate points are calculated,® it only remains to 
compare their minimal distances to the existing points p* G P and then, take the one 
whose distance is maximal; in fact, the resulting set is hnite and reasonably small. 
This is the CIE76 Voronoi method of solving Problem 1 with respect to the Euclidean 
distance. In Table 2 we list the hve new colors that we obtained, one by one, for the 
Moscow metro by using this method. 

It is readily seen that the Voronoi method would provide the global optimal solution 
to maximin Problem 1 — if there were a perceptually uniform color space. This is not 
yet the case nowadays, therefore let us study how better extimates can be obtained. 

Because the CIEDE2000 distance is much more accurate than the Euclidean distance, 
we can improve our method by selecting the optimal point from the candidate points 
by comparing the APqo distances instead of the Euclidean distances XEah- This is the 
combined CIE76-CIEDE2000 Voronoi method. Of course, this method does not solve 

^Metro lines are built one after another, so choosing new colors one by one is quite logical. However, 
sometimes we know in advance that several new lines are being built. 

®This can be done by elementary linear algebra, that is, by calculating bisectors and intersections 
of planes. In Figure 4, a 2-dimensional example of a Voronoi diagram can be found. 
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Figure 4. A 2-dimensional Voronoi diagram. Note that all edges of the 
Voronoi polygons are parts of the bisectors of lines connecting two points 
in the diagram. 


Color 

min AEab 

min AEqo 

L 

a 

b 

• Bright green 

114,4 

22,0 

83 

-138 

91 

• Blue 

87,0 

19,8 

33 

80 

-109 

Sea green 

67,7 

13,1 

85 

-106 

31 

Cyan 

65,3 

21,9 

87 

-78 

-21 

• Bright pink 

56,1 

14,4 

60 

100 

-64 


Table 2. Colors we hnd one by one for the first five new lines in the 
Moscow metro by using the CIE76 Voronoi method, that is, by solving 
Problem 1 with respect to the Euclidean distance. The minimal distances 
which we indicate here take in account the points that came earlier into 
the table. Note that the minimal AEqo distance is then not necessarily 
decreasing. 


Problem 1 with respect to AEqo alone, because the candidate points are found by using 
AEab- However it is a fast method to obtain colors by using the CIEDE2000 distance. 
The five new colors obtained, one by one, for the Moscow metro by this method are 
given in Table 3. We see that this gives us much better solutions with respect to the 
CIEDE2000 distance than the Voronoi CIE76 method: indeed, the AEqo distance drops 
at a much slower rate. 


5. The simplex method 

We can also solve Problem 2 with respect to the Euclidean distance by using the simplex 
method. Let us consider here the case when the number of colors m is two, which gives 
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Color 

min AEqq 

L 

a 

b 

Cyan 

23,9 

87 

-78 

-21 

• Bright green 

22,0 

83 

-138 

91 

• Blue 

19,8 

33 

80 

-109 

Pale lavender 

19,7 

86 

23 

-22 

• Rose pink 

19,8 

78 

59 

12 


Table 3. Colors we find one by one for the first five new lines using the 
combined CIE76-CIEDE2000 Voronoi method. Note that the distances 
here are the CIEDE2000 distances, and they drop at a much slower rate 
than in Table 2.® 


us the following function to maximize, 

M{x,y) = min{||a: -Pill, \\y - Pi\\, ||x - i/||}. 

Pi&P 

x,y&T 

We have (x, p) G E x E C M®. Let us rewrite our problem: hrst, put one extra parameter 
xo and the two vectors x,y together in x = {xo,x,y) G and now maximize the 
function Mq{x) = xq over all x G M x E x E satisfying the constraints 


1 X — Pi| — Xo > 0 

for each pt G P, 

(la) 

lb-Pill -a^o > 0 

for each pj ^ P, 

(lb) 

||x - y\\ - Xo > 0. 


(Ic) 


Of course, more constraints on (x, y) are used to delimit the set E x E. Now we use the 
simplex method to maximize our function. 

The gamut’s boundary gives the linear constraints, and the nonlinear constraint 
functions are equations (1). Because the simplex method gives us only a local optimum, 
we apply it for a million random initial points (xo,x, p) G (with xq G [0,20] and 
{x,y) G E X E). These are local optima with respect to the AEab distance; let us now 
choose the one with maximal AEqo distance. This way, we obtain the two colors 

Bright aqua (84, -65, -12) and Baby pink (86, 35, 7). 

Their combined minimal CIEDE2000 distance to the old colors and each other is 22,7. 
This is more than the minimal distance which was obtained after two steps of the 
one by one CIE76-CIEDE2000 Voronoi method. So this gives us a better result for 
solving Problem 2. However, this method is quite time-consuming (as m increases, the 
simplex method will need even more time), so the much faster Voronoi method is a 
good alternative. 


6. Conclusion and discussion 

We conclude that the Voronoi diagram method (and its improvement by using both 
the CIE76 and CIEDE2000 distance functions) for adding new colors one by one gives 

®The minimal CIEDE2000 distances do not always decrease each step; this is because the candidate 
points are selected by using the Euclidean distance, which does not necessarily give an optimum with 
respect to the CIEDE2000 distance. Therefore it is possible that, after adding a new color, one of the 
new candidate points is even a better solution than the point added before. 
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solutions which are comparable with the ones produced by the nonsmooth optimization 
method for hnding two or more new colors at once. This is true for the map of the 
Moscow metro under study. Other solution schemes are of course possible, e.g. the 
Monte-Carlo method: test points Xj G T for the new color(s) are generated within 
the gamnt at random, and the new color - or a tuple of new colors - is the maximin 
optimum over a sufficiently large number of such attempts. The ballistic method is 
realized as follows: first, suppose that each old color Pi G T and every new point Xj in 
the sought-for tuple carries a positive electric charge (so that they repulse by Coulomb’s 
law). To prevent an escape of new points from the gamut T across its boundary, put a 
negligibly small but positive charge at the boundary’s closest point to Xj for every j. 
By starting from random initial data a;j(0) G T, the conhgurations {xj{t)} of charges 
evolve inside the polyhedron T which, by assuption, is hlled with a viscous medium so 
that the new points’ motion slows down; sooner or later they are captured near the 
local minima of the potential (see [1] for details on both methods). 

Finding several new colors is still not all that one needs in order to paint the new 
lines on a metro map. Clearly, it is the conhguration of transfers from the old to new 
lines that must further be taken into acconnt — to distribute the new colors between 
the new lines in snch a way that at every transfer station, the colors of intersecting lines 
are maximally different. 

Breaking news: an on-line poll. The 16.8 km long Kozhnkhovskaya line with eight 
stations on it is scheduled to appear in the Moscow undergronnd system in 2015-16, 
officials say. The new line will connect the (sonth-)east with the city centre. 

Two colors, to freely choose from, were offered to pre-registered “active citizens” 
dnring an interactive poll, held from 20 October till 10 November 2014. Namely, black 
and pink were claimed to remain the only colors from the gamut that were not yet used 
to paint lines on the Moscow metro map (see Fignre 1). As many as 307,350 active 
citizens took part in the on-line poll, media report; 71% voted pink and 18% preferred 
black, whereas 7% trusted the choice to experts (such as designers, ethnographers, or 
psychologists). Finally, 4% of the respondents either made their own suggestions of the 
new color: tnrqnoise, olive, coral, etc. - or even proposed to make that line on the map 
speckled (see [10, 11]). 

Permitting a large group of people (here, 307,350 pre-registered active citizens) to 
choose the new color(s) from a given list of pre-selected variants (here, two: black vs 
pink) might not always be the optimal strategy to resolve an important issne, painting 
new line(s) on a metro map. Perhaps, the decision should be trusted to the Protector of 
tunnel builders and Best Friend of all cartographers; then, the hnal choice is sometimes 
delegated to a child. Alternatively, one could just ask a mathematician. 
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